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ABSTRACT 

We show that the 5-GHz radio flux of transient ballistic jets in black hole binaries corre- 
lates with the dimensionless black hole spin parameter a* estimated via the continuum-fitting 
method. The data suggest that jet power scales either as the square of a* or the square of the 
angular velocity of the horizon £Ih ■ This is the first direct evidence that jets may be powered 
by black hole spin energy. The observed correlation validates the continuum-fitting method of 
measuring spin. In addition, for those black holes that have well-sampled radio observations 
of ballistic jets, the correlation may be used to obtain rough estimates of their spins. 
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1 INTRODUCTION 

Accreting black holes, both supermassive and stellar m ass, are 
commonly observed to produce p owerful relativistic jets dZensus! 
1 19971 : iMirabel & Rodriguez! 1 19991) . Although there now exists a 
wealth of data and many detailed models, the mechanism that pow- 
ers these jets remains a mystery. 

The popular idea th at jets are power ed by the black hole (BH) 
goes back to the work of lPenrosej (1 19691), who s howed that a spin- 
ning BH has free energy. iBlandford & Znaiekl d 19771) proposed a 
plausible mechanism whereby this free energy could be used to 
power an astrophysical jet. They suggested that magnetic fields in 
the vicinity of an accreting BH would be twisted as a result of the 
dragging of space-time by the rotating BH. The twisted field lines 
will carry away energy from the BH, producing an electromagnetic 
jet. The broad outlines of this model have been co nfirmed in nu- 
merical simulations (e.g. lTch ekhovsko y et alj|201ll) . 

While a connection between BH spin energy and relativistic 
jets is theoretically very appealing, there has been no direct obser- 
vational evidence for such a link. This is because, until recently, 
there was no BH with a believable measurement of the dimension- 
less spin parameter et» = cJ/GM 2 , where M and J are the mass 
and angular momentum of the BH. The situation has now changed. 
Methods are now avai lable — one in pa rt icular, the continuum- 
fitting (CF) method 1 ! dZhang et all Il997l: iGierlinski et alj I2OO1I : 
IShafee et alj|2006l : iDavis et all 120061 : iMcClintock et alj|2006|) — 
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1 A second method, based on modeling the relativistically broadened X-ray 
iron Ka line, is not considered in this Letter (see Section 4 for a discussion). 



that have enabled us to make plausibly reliable measurements of 
a, for several stellar mass BHs. With this sample of spin measure- 
ments, we are now in a position to test whether jet power is related 
to BH spin. Such a test is the goal of this Letter. 

In Section 2, we describe our sample of stellar mass BHs and 
collect together the relevant observational data on BH spins and 
jet power. In Section 3, we plot radio power against BH spin and 
demonstrate that there is a significant correlation between the two 
quantities. We summarize and discuss in Section 4. 



2 THE DATA 

2.1 BH sample and spin estimates 

The CF method (see lMcClintock et alj|201 ll for a brief review) fits 
the X-ray continuum spectrum of an accreting stellar mass BH us- 
ing th e classic relativistic thin-disc model of iNovikov & Thornd 
dl973l) . The spectral fit gives an estimate of the radius of the in- 
ner edge of the accretion disc. The BH spin parameter a» is then 
obtained by assuming that the disc edge is located at the innermost 
stable circular orbit (ISCO) of the Kerr metric. The CF method has 
been developed in detail over the last several years and has been 
shown to produce consistent results when multi ple independent ob- 
serva tions of the same source are available (e.g. lsteiner et alj2009l 
120101) . In addition, numerical simulations have provided support 
for a crucial assumption of the model, namely that the disc inner 
edge is close to the ISCO dShafee et al.ll2008l : IPenna et alj|2010l : 
iKulkarni et al.11201 lMNoble et al.11201 ll) 

The spins of the BH primaries in ni ne BH binaries (B HBs) 
have been measured using the CF method dMcClintock etal.l201ll : 
lGouetalj|20TTh . Five of these BHBs, namely, A0620-00, XTE 
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Table 1. Parameters of transient BHBs in the sample 



BH Binary a, M (Mq) D (kpc) i (deg) (&)max,5GHz (Jy) S (7 = 2) (Jy) References 



A0620-00 


0.12 ±0.19 


6.61 ±0.25 


1.06 ±0.12 


51.0 ±0.9 


0.203 


0.145 


1,6,7 


XTEJ1550-564 


0.34 ±0.24 


9.10 ±0.61 


4.38 ±0.50 


74.7 ±3.8 


0.265 


0.859 


2, 6,8 


GRO J 1655^10 


0.7 ±0.1 


6.30 ±0.27 


3.2 ±0.5 


70.2 ± 1.9 


2.42 


7.74 


3, 4, 6, 9, 10 


GRS 1915+105 


0.975 ±0.025 


14.0 ±4.4 


11.0 ± 1.0 


66.0 ±2.0 


0.912 


2.04 


5,6, 11, 12 


4U 1543^17 


0.8 ±0.1 


9.4 ± 1.0 


7.5 ± 1.0 


20.7 ±1.5 


> 1.16 x 10~ 2 


> 4.31 x 10~ 4 


3, 6, 13 



Re ferences: (nlGou et aljl2010l): (2)ISteiner et alJfeOl ll) : (3) |Shafee et al.N2006l), (4)|Pavis et alJbOOrj): (5)lMcCfintock eTafl 120061): (6)l6zel et al.N201(j) ; 
(7)lKuulkers et all h999H; (8)lHannikainen et alj d2009l) : (9) iHiellming & Rupenl Il995h ; (lO) lHannikainen et alj feOOOl) ; (1 P lRodriguez et al.Ul995l) : (12) 
lFenderetalJil999l) ; (13) |Park et aljj2004) . 



J1550-564, GRO J 1655^0, GRS 1915+105, 4 U 1543-47, are 
transient systems l lRemillard & McClintockll2006t) . These five sys- 
tems have low-mass secondaries and undergo mass transfer via 
Roche lobe overflow. They are of primary interest to us because 
during outburst, as they approach the Eddington limit, they produce 
ballistic jets (Section \2.2\ . The measured BH spin values a* and 
masses M, along with distances D and binary inclination angles i, 
are listed in Table 1 . In the case of A0620-00 and XTE J 1 550-564, 
the error estimat e s on the spins are t aken from the original papers 
dGou et al1l2010l : ISteiner et alj|201 lh . The oth er three spins were 
measured in the early days of the CF met hod l lShafee et alj|2006t 
iDavis et al. I l2006l : iMcClintock etal]|2006h . and we have doubled 
the published error estimates. 

An additional fou r stellar ma ss BHs have spin estim ates: LMC 
X-3 jDavis et al.l2006h . M33 X- 7 dLiu et al.l2008l . l2010l) , LMC X-l 
dGou et al.ll2009h and Cyg X-l dGou et alfcoill) . These are persis- 
tent BHBs dRemillard & McClintockll2006|)~ which have high-mass 
companion stars and undergo mass transfer via winds. Also, they 
do not show the kind of transient behavior seen in the previous five 
objects and are generally understood to belong to a different class. 
We ignore them in this study. 
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2.2 Jet radio power 

iFender et alj d2004l) identified a number of systematic properties 
in the radio emission of BHB jets. They showed that there are 
two kinds of jets associated with specific spectral states of the 
X-ray source. The first type of jet is observed in the hard spec- 
tral state as a steady outflow. This jet is observable only out to 
a few tens of au and is apparently not very relativistic. The sec- 
ond and far more dramatic jet, which is central to this Letter, is 
launched when a BHB with a low-m ass secondary undergoes a 
transient outburst dFender et alll2004l) . This powerful transient jet 
usually appears near (or soon after) the time of outburst maxi- 
mum, as the source switches from its inital hard state to a soft state 
via the 'steep power-law' (SPL) state, a violently-variable state 
characterize d by both strong thermal and po wer-law components 
of emission dRemillard & McClintockll2006l) . Transient jets man- 
ifest themselves as blobs of radio (and occasionally X-ray) emit- 
ting plasma that move ballistically outward at relativistic speeds 
(7j ct > 2). Because these pc-scale jets resemble the kpc-scale jets 
seen in quasars, BHBs that produce them are called microquasars 
dMirabel & Rodrigued[l999h . 

Ballistic jet ejection occurs at a very specific stage during the 
spectral evolution of a given system dFender et al1l2004h . As most 
clearly demonstrated for the prototype microquasar GRS 1915+105 
dFender & Bellonill2004h . this stage appears to correspond to the 
inward-moving inner edge of the accretion disc reaching the ISCO, 



Figure 1. Plot of the maximum observed radio power (vSv)mz.x as a func- 
tion of frequency v for transient ballistic jet outbursts in four BHBs. Two 
separate outbursts are shown for GRS 1915+105. Best-fit lines (two sepa- 
rate ones in the case of GRS 1915+105) are indicated, except in the case of 
A0620-00 where the line slope is fixed at 0.6 (or a = —0.4). 

which results in a shock or some other violent event that launches 
the large-scale relativistic jet. In this scenario, it appears reason- 
able that the jet is launched within a few gravitational radii and 
hence plausible that the spin energy of the BH could power the 
jet. In contrast, the steady jet in the hard stat e is thought to origi - 
nate much further out at ~ 10 - 100 GM/c 2 dMarkoff et aHl2005t) 
where the effects of spin are relatively weak. Another virtue of the 
ballistic jets for our purposes is that they occur at a sharply defined 
luminosity (i.e. near Eddington) compared to the hard state steady 
jets, which occur over a wide range of luminosity. Ballistic jets are 
thus better 'standard candles'. For these reasons, in this Letter we 
restrict our attention to ballistic jets from transient low-mass BHBs. 

A typical ballistic jet blob is initially optically thick and has 
a low radio power. As the blob moves out and expands, the larger 
surface area causes its radio power to increase. This continues un- 
til the blob becomes optically thin, after which the flux declines 
rapidly. The overall be haviour is generally consisten t with an ex- 
panding conical jet (e.g. IHiellming & Johnstoniri988l) . 

Fig. shows the peak radio flux (5 , I y) max versus v observed 
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at different radio frequencies v for four of the five transient BHBs 
in our sample. The radio light curves of these four systems were 
monitored with good time resolution, allowing us to obtain rea- 
sonably accurate estimates of the peak fluxes. The top left panel 
shows data for two separate outbursts of GRS 1915+105 (the solid 
an d open circles corresp ond respectively to the outbursts studied 
bv lRodriguez et al J 19951 and lFender et all 19991) . The two lines are 
fits to the respective data and have a slope of 0.59; writing the spec- 
trum as S v oc u a , the fit correspon ds to a = —0.41. The top right 
panel combines the obser vations o flHiellming & Rupenl fl99a) and 
lHannikainen et alj < l2000h during an outburst of GRO J 1655^10. 
The best-fit line corresponds to a = —0.663 The lower two 
panels show data for XTE J 1550-564 dHannikainen et al.l 120091 . 
a = -0.18) and A0620-00 l lKuulkers et alJI 19991) . FoTthe latter 
source, we do not have enough data points to determine the slope; 
the line in the plot corresponds to a = —0.4, the average spectral 
index of the other three BHBs. In order to enable a fair comparison 
of the different objects, we use the fitted lines in the four panels to 
estimate the peak fluxes (S^max at a standard frequency of 5 GHz. 
These 5-GHz peak flux values are listed in Table 1. 

While each of the above four objects was densely observed 
in radio during one or more transient outbursts, 4U 1543-47 was 
unfortunately not monitored well at radio frequencies during any 
of its several outbursts. The only radio data we know of when the 
source was brigh t are those for the 2002 outburst summarized in 
IPark et alj <2004 . The strongest radio flux was 0.022 Jy at 1 .02675 
GHz. Assuming a = —0.4, this gives a flux of 0.01 16 Jy at 5 GHz 
(or only 0.00043 Jy if one corrects for beaming with 7j ct = 2). We 
list this result separately in Table 1 and plot it as a lower limit in 
Figs 2 and 3 because of the sparse radio coverage. In addition, there 
was an anomaly in the 2002 X-ray outburst of this source. 

The anomalous beh aviour of 4U 1543^-7 is apparen t by an 
inspection of figs 4-9 in iRemillard & McClintockl 120061) . which 
summarize in detail the behaviour of six BH transients scrutinized 
by RXTE. In panel b of these figures, which displays light curves 
of the PCA model flux coded by X-ray state, one sees that only 
4U 1543-47 failed to enter the SPL state (green triangles) near the 
peak of its outburst, i.e. at the time of the radio coverage reported 
by Park et al. Rather, it remained locked in the thermal state (red 
crosses) after its rise out of the hard state. This behaviour contrasts 
sharply with the behaviour of the other five transients which dis- 
played the strongly-Comptonized SPL state during both the late 
phase of their rise to maximum and during their early decay phase. 
Thus, because of (1) the sparse radio coverage of 4U 1543-47, and 
(2) the failur e of the source to transition out of the jet-quenched 
thermal state dGallo et alJPiool) to the SPL state (which is closely 
associated with the launching of ballistic jets), we treat the maxi- 
mum observed flux of 0.022 Jy as a lower limit. Fin ally, in sharp 
contra st to our finding, we note that figs 5 and 6 in lFender et alj 
d2004l) indicate a very high jet power for 4U 1543-47. We are un- 
sure how they arrived at their result, but suspect it was based on 
infrared data and their equipartition model (see Section 4). If so, an 
extension of the present work to infrared data might be worthwhile. 

To measure jet power, we scale the 5-GHz peak flux of each 
BHB by the square of the distance to the source D. We also divide 
by the BH mass M since we expect the power to be proportional 
to M (this scaling is not important since the range of masses is 



2 In the case of GRO 11655^10, the 22-GHz observations did not cover the 
peak of the light curve. Hence this point is shown as a lower limit. Similarly, 
in A0620-00, the peak was not observed at 0.962 and 1.14 GHz. 
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Figure 2. Plot of the jet power Pj ot as estimated from the maximum radio 
flux of ballistic jets (equation 1) vs the measured spin parameter of the BH 
a* for the transient BHBs in our sample. Solid circles correspond to the first 
four objects listed in Table 1, which have high quality radio data, and the 
open circle corresponds to 4U 1543^7, which has only a lower limit on the 
jet power. The dashed line corresponds to Pj ct oc a 2 , the theoretical scaling 
derived bv lBland ford & Znaiekl 119771) . The data suggest that ballistic jets 
derive their power from the spin of the central BH. 



only a factor ~ 2). We thus obtain from the radio observations the 
following quantity, which we treat as a proxy for the jet power: 



(1) 



It is hard to assess the uncertainty in the estimated values of Pjct- 
There is some uncertainty in the values of D and M, but these are 
not large. Potentially more serious, the radio flux may not track 
jet power accurately. For instance, the properties of the ISM in the 
vicinity of the BHB may play a role and are likely to vary from 
one object to another. Also, the energy released in these roughly 
Eddington-limited events will vary (e.g. see GRS 1915+105 in 
Fig. 1). Below, we arbitrarily assume that the uncertainty in Pj c t 
is 0.3 in the log, i.e. a factor of 2 each way. 



3 JET POWER VS BH SPIN 

Fig. 2 shows jet power Pj c t plotted against BH spin parameter a* 
for the five transient BHBs in our sample. The data are taken from 
Table 1 . The dashed line has a slo pe of 2, motivated by t he the oreti- 
cal scaling, Pj ct oc a„, derived bv lBlandford & Znaiekl dl977l) . The 
data points agree remarkably we ll with this theoretical prediction. 

iBlandford & Znaiekl d 19771) assum ed a slowly spinning BH: 
a* <C 1. iTchekhovskov et alj Hoioh obtained a more accurate 
theoretical scaling which works up to spins fairly close to unity: 



Pjct oc Q. H , where Qh is the angular frequency of the BH hori- 
zon, fl H = a,(c 3 /2GM)/(l + a/1 - al). Fig. 3 shows a plot of 
Pjet vs SI if , with the dashed line corresponding to a slope of 2. The 
agreement is again very good. 
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Figure 3. Similar to Fig. 2, but showing the angular velocity of the BH 
horizon along the absc issa. The dashed line corresponds to Pj et oc 
iTchekhovskov et alj20Toh . 



We need to consider one additional effect: relat ivistic beam- 
ing. A ssuming a typical jet Lorentz factor 7j ot = 2 dFender et al.l 
120041) and using the inclination angles i given in Table 1, 
we have corrected the values of (S v ) ma ,x,5 GHz- The beaming- 
corrected radio fluxes So ( computed using the relations given in 
Irvlirabel & Rodriguez! 1 19991 with the values of a given in Section 
2.2) are listed in Table 1. The inferred jet power of A0620-00 (the 
left-most point in Figs 2 and 3) decreases by a small factor, whereas 
the other three jet powers increase by a larger factor. As a result, 
the spread in Pj ct among the four objects becomes ~ 3.0 orders of 
magnitude, compared to ~ 2.4 orders in Figs 2, 3. Thus, allowing 
for beaming enhances the range of Pj ct in the sources and solidifies 
further the connection between jet power and BH spin. 



4 DISCUSSION 

Since the correlations shown in Figs 2 and 3 are based on only four 
objects, one wonders whether we are merely seeing chance align- 
ment of intrinsically uncorrelated data. The chief argument against 
this hypothesis is that a* varies over the full allowable range of 
prograde spins, Qh varies by more than a factor of 10, and Pj ot 
varies by 2.4 orders of magnitude (or 3 orders of magnitude if one 
corrects for beaming assuming 7j e t = 2). Also, the plotted points 
differ from one another by several standard deviations, which is not 
statistically likely. Therefore we conclude that the power of ballis- 
tic jets is most likely correlated with the spin of the accreting BH. 

At the same time, and as a corollary, the strong apparent corre- 
lation validates the CF method of measuring spin. The CF method 
is based on a number of ass umptions, most of whic h have been 
independently validated (see iMcClintock et al .11201 it) . The results 
presented here provide yet another validation. Caveats to the above 
conclusions include the small size of the sample, insufficient data 
on one object (4U 1543-47), and uncertainties in how well jet 
power and radio luminosity track each other. 



The existence of a correlation between jet power and BH spin 
does not necessarily mean that the energy source for the jet is BH 
spin. The power could possibly be supplied by the accretion disc 
dGhosh & Abramowiczlll997l : lLivio et al.ll 1 9991) . Since the binding 
energy of a particle at the ISCO increases with increasing a*, the 
disc power increases with BH spin and this might cause the ob- 
served correlation. However, we note that the radiative luminosity 
of a thin accretion disc varies by only a modest factor with BH 
spin; the radiative efficiency rj = 0.061 for a* = 0.12 (the spin 
of A0620-00) and 77 = 0.23 for a* = 0.98 (the spin of GRS 
1915+105). If jet power scales similarly, and if radio luminosity is 
roughly proportional to jet power, we expect no more than a factor 
of 4 variation in the radio powers in our sample. Instead, we see a 
factor ~ 250 — 1000. Moreover, the observed spread is rather close 
to what is expected theoretically if jets are powered by BH spin. 
The evidence thus suggests that ballistic jets are powered directly 
by the spin energy of the accreting BH. 

Based on the above arguments, we view our results as a confir- 
mation of the Penrose-Blandford-Znajek mechanism of powering 
relativistic jets by BH spin energy. Theoretically, this mechanism 
depends on both the BH spin and the magnetic field strength at 
the horizon. The latter is believed to d epend on the mass accretion 
rate M (e.g. ITchekhovskov et alj201ll) . Since ballistic jets are seen 
during a very specific phase of the evolution of a transient BHB, it 
is reasonable to assume that M (normalized by the BH mass) is 
roughly the same in different objects when they exhibit ballistic 
jets, or indeed in different ballistic jet episodes in the same object. 
This allows us to eliminate M from our analysis and to treat bal- 
listic jets as 'standard candles', thereby making the comparisons 
shown in Figs 2 and 3 meaningful. 

In addition to the CF method, a second method based on fit- 
ting the profile of the relati vistically broadened Fe Kg line has bee n 
used to estimate BH spins jRevnolds & Nowakl2003l ; lMillerl2007l) . 
In this Letter, we opt to use only CF spin data for two reasons: (1) 
The Fe-line models are complex and therefore relatively less re- 
liable. CF spins are obtained rather simply by modeling a domi- 
nant thermal disc component, while Fe-line spins require modeling 
the thermal disc plus a Compton component plus a disc reflection 
component, which includes the Fe Kq line. The Fe-line method 
furthermore requires characterizing a luminous corona of unknown 
geometry. (2) For several well-studied systems, the Fe-line method 
has generated widely inconsistent values of the spin parameter or 
show n to be strongly m odel dependent (e.g. fo r Cyg X-l, see Sec . 

7.1 in lGouetalJboill ; for MCG-6-30-15, see iMiller et al]|2009h . 
The CF method, on the other hand, gives consistent results for mul- 
tiple and independent observations of individual sources. For ex- 
ample, for the BHBs listed in Table 1 (excluding A0620-00), con- 
sistent results were obtained for « 50 RXTE spectra (XTE J 1550- 
564); 2ASCA and 31 RXTE spectra (GRO J1655-40); 1 ASCA and 
5 RXTE spectra (GRS 1915+105); and 34 RXTE spectra (4U 1543- 
47). The standout example is LMC X- 3 with 411 spectra collected 
by eight X-ray missions over 26 years dSteiner et alfeOldt) . 

After considering s eparat ely both ballistic and hard state 
steady iets. lFender et alj d2010h find no evidence for a correlation 
between jet power and BH spin. We have already given in Section 

2.2 a plausible reason for the absence of evidence in the case of 
the steady jets. We now focus on rationalizing the very different re- 
sults obtained by Fender et al. (presented in their Section 2.2.1 and 
fig. 6) and ourselves for the ballistic jets. Our data sample (Table 1) 
is identical to their comparable sample (see the right panel of their 
fig. 6). The only significant difference in dat a selection is that fo r 
GRO J 1550-564 we use the new spin value of lSteiner et al. 
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a, - 0.34 ± 0.24, while they used the earlier lDavis et alj d2006h 
limit of a* < 0.8. 

The substantia l difference between our results and those of 
iFender et alj feOld) is, in the end, determined by the choice of 
the quantity used to represent jet power. We simply use the max- 
imum observed flux density at 5 GHz expressed as a luminosity. 
Fender et al. compute jet power from the peak radio luminosity 
and the rise time of some particular synchrotron event. The authors 
clearly state that their approach 'is useful to provide lower lim- 
its on, and order-of-magnitude estimates of, jet power but is very 
susceptible to errors resulting from poor sampling of events, un- 
certainties in Doppler boosting, assumptions about equipartition, 
etc' Their es timates of jet lumino sity for three sources are given 
in Table 1 of IFender et al.1 d2004l) , but it is not clear how the lu- 
minosities of A0620-00 and 4U 1543-47 were estimated. The au- 
thors further adopt a formula relating jet power to X-ray luminosity, 
logiol/jct = c + 0.5(logiol/ x — 34), and estimate the normaliza- 
tion constant c in the preceding formula, which they treat as their 
proxy for jet power. In short, their proxy for jet power is heavily 
model dependent and ours is model independent. 

The correlation shown in Fig. 2 can be used to obtain rough 
estimates of spin for any transient BH that has undergone a ma- 
jor outburst cycle and that has been closely monitored at radio 
wavelengths. For instanc e, radio observati ons of Nova Muscae 
1991 (GRS 1124-68) by iBall et al] ( 1 19951) suggest a maximum 
5-GHz radio flux ~ 0.2 Jy. Assum ing a distance D ~ 6 kpc 
and a typical BH mass M m 8M Q \6zel et al.|[2oToh . we obtain 
log[Pj c t] ~ 0.65±0.3. Fig. 2 then gives a, « 0.3 — 0.6 . In the case 
of GX 339-4, the brightest X-ray and radio outburst jGallo et al] 
l2004h had a maximum 5-GHz flux of 0.055 Jy. Taking D « 9 kpc, 
M « 8M dOzel et al.|[2oToh . we find log[P jct ] « 0.4 ± 0.3 and 
a, ~ 0.2 — 0.5. The latter esti mate is consistent with the str ict up- 
per limit a, < 0.9 derived bv lKolehmainen & Done! feOlCh using 
the CF method with conservative assumptions. 

These examples illustrate the importance of obtaining good ra- 
dio coverage for all future transient BHBs, including especially the 
recurrent system 4U 1543-47. Those systems that have CF-based 
spin measurements will flesh out the correlations plotted in Figs 2 
and 3. For the many other BH transients that lack a sufficiently- 
bright optical counterpart and are therefore out of reach of the CF 
method, the radio data can either be used as a check on Fe-line 
measurements of spin or serve as our only estimate of spin. 
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